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Abstract 



We calculate the tunneling conductance spectra of a ferromagnetic 

metal/insulator /superconductor using the Blonder-Tinkham-Klapwijk (BTK) 

formulation. Two possible states for the superconductor are considered with 

the time reversal symmetry (T) broken, i.e., d^2_y2 +is, or d^2_y2 +idxy In 

both cases the tunneling conductance within the gap is suppressed with the 

increase of the exchange interaction due to the suppression of the Andreev 

reflection. In the {dx2_y2 + zs)-wave case the peaks that exist when the fer- 

romagnet is a normal metal in the amplitude of the s-wave component due 

to the bound state formation are reduced symmetrically, with the increase of 

the exchange field, while in the {dx2_y2 + idxy)-wave case the residual density 

of states within the gap develops a dip around E = with the increase of 

the exchange field. These results would be useful to discriminate between 

T-broken pairing states near the surface in high-Tc superconductors. 
74.20. z, 74.50.+r, 74.80.Fp 
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I. INTRODUCTION 



It is generally accepted that the pairing state of high-T^ cuprate superconductors is of 
dominant {dx2_y2)-wa.ve symmetry. However, it is possible for a secondary component to 
be induced wherever the dominant order parameter varies spatially, for example, in regions 
close to surfaces |Q. The local pairing symmetry can be detected by tunneling measurements 
where Andreev reflection take place [0,^. In the Andreev reflection process an electron in- 
cident, in the barrier with an energy below the superconducting gap can not drain off into 
the superconductor. It is instead reflected as a hole and a Cooper pair is transferred into 
the superconductor. In anisotropic high-Tc superconductors the transmitted quasiparti- 
cles experience different sign of the pair potential. This results in the formation of bound 
states close to the surface @], which are detected as zero energy peaks in the tunneling 
spectra as an indication of {da;2_y2)-waMe pairing symmetry ^j. The splitting of the zero 
energy conductance peak at low temperatures observed experimentally [|^,|^ is a signature 
of a two component pairing state with the time reversal symmetry broken, and is consis- 
tent with both the ci^z-ya + is, or d^2_y2 + id^y pairing state. The tunneling spectra of a 
normal metal/insulator/superconductor junctions, for the above pairing states has already 
been calculated, by extending the Blonder-Tinkham-Klapwijk (BTK) theory to include the 
anisotropy of the pair potential Based on the bound state formation, simple arguments 
have been derived to explain the subgap conductance in each case and to discriminate be- 
tween the two pairing symmetries. 

Moreover the tunneling spectra of ferromagnet /superconductor junctions has been clar- 
ified both for singlet fl-[r2|, and triplet |13| pairing states where the important parameter 
is the exchange field. Also the changes of the critical temperature of the superconductor in 



s-wave-superconductor- ferromagnet multilayers have been studied [|T4[. In this paper we use 
the scattering approach to study the charge transport through ferromagnet /superconductor 
with time reversal T-breaking pairing states junctions. For the superconductor we assume 
two pairing states i.e. dx2-y2 + is, or dx2-y2 +idxy, with a complex order parameter breaking 
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the time reversal symmetry. In both cases the tunnehng conductance within the gap is sup- 
pressed with the increase of the exchange interaction, due to the suppression of the Andreev 
reflection. In the {d^2_y2 +zs)-wave case the peaks that exist when the ferromagnet is a nor- 
mal metal in the amplitude of the s-wave component due to the bound state formation are 
reduced symmetrically. In the ((i^2_y2 + idxy)-wave case the tunneling conductance within 
the gap develops a dip around E = with the increase of the exchange field. The tunneling 
conductance characteristics for T-broken pairing states are discussed extensively in Sec. III. 

The evolution of the Andreev- and normal-reflection amplitude with the exchange field 
presented in Sec. IV explains the suppression of the tunneling conductance with the exchange 
field. The Andreev reflection amplitude decays to zero at a critical value of the exchange 
field Xc that depends only on the angle 9 of the quasiparticle trajectory. A bound state 
at 6 contributes to the Andreev reflection and to the conductance for a given x only when 

X *\ X Q • 

The magnetic field for finite exchange interaction induces an imbalance in the peak 
heights, of the spectra for energy positive and negative as seen in Sec. V. 

Also the c-axis tunneling spectra presented in the last section is similar for T broken 
pairing states, since the pairing potential along the c-axis does not change sign on the Fermi 
surface. These results would be helpful to discriminate between time-reversal symmetry 
broken surface states in high superconductors. 



II. THEORY OF TUNNELING EFFECT 

The motion of quasiparticles in inhomogeneous superconductors is described by the Bo- 
goliubov de Gennes (BdG) equations. The motion of electrons inside the ferromagnet is 
described within the Stoner model by an effective single particle Hamiltonian with an ex- 
change interaction. The BdG equations read [|T^] 

{Heir)- pU{x))u{r) + Sdr'A{s,x)v{r') = Eu{r) 
J dr'A*{s, x)u{r') - {HUr) + pU{x))v{r) = Ev{r) 
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The single-particle Hamiltonian is given by Heir) = -h'^ Vr /2"^e + V{r) - Ep, E is the 
energy measured from the Fermi energy Ep. U{r) is the exchange potential, p is 1(— 1) 
for up(down)-spins. A{s,x) is the pair potential, after a transformation from the position 
coordinates r,r' to the center of mass coordinate x = (r + r')/2 and the relative vector 
s = r — r'. After Fourier transformation the pair potential depends on the related wave 
vector k and x. In the weak coupling limit k is fixed on the Fermi surface = kp), and 



only its direction 9 is variable. After applying the quasi-classical approximation, i.e., |jT5| 



^ u(r) ^ 



-ikr 



( uir) ^ 



(2) 



so that the fast oscillating part, of the wave function is divided out, the BdG equations are 
reduced to the Andreev equations H] 



(3) 



Eu{r) = —ivpk • \/u{r) + A{9,r)v{r) 
Ev{r) = ivpk ■ \/v{r) + A*{9,r)u{r) 

where the quantities u{r) and v{r) are electronlike and holelike quasiparticles and vp is the 
Fermi velocity. 

We consider the ferromagnet /insulator /superconductor junction shown in Fig. 1. We 
choose the y direction to be parallel to the interface, and the x direction to be normal to 
the interface. The insulator is modeled by a 5 function, located at x = 0, of the form V6{x). 
The temperature is fixed to K. We take both the pair potential and the exchange energy 
as a step function i.e. A{6,r) = Q{x)A{6), U{r) = Q{—x)U. For the geometry shown in 
Fig. Eqs. 1^ take the form 



(4) 



Eu{x) = -ivFkp^£u{x) + A{e)v{x) 
Ev{x) = ivpkpx^v{x) + A*{6)u{x) 

When a beam of electrons is incident from the ferromagnet to the insulator, with an 
angle 9, the general solution of Eqs. (Q), is the two component wave function that for 
a; < is written as 



4 



^1 



^iq l[t]X COS 9 A _j_ Jj 



(A 



(5) 



where a^[|], fe-f^j], are the amphtudes for Andreev and normal reflection for spin- up (-down) 
quasiparticles, and = 



^{Ep ± U) is the wave vector of quasiparticles in the ferro- 
magnet for up (down-)-spin. The wave vector of the electronlike, holelike quasiparticles is 
approximated by ks = \J ^™f'^ ■ Since the translational symmetry holds in the y-axis direc- 
tion, the momenta parallel to the interface is conserved, i.e., sm9 = sin^^ = ksSinOg. 
Note that 6 is different than 6a since the retroreflection of the Andreev reflection is bro- 
ken. Using the matching conditions of the wave function at x = 0, \E'/(0) = \E'//(0) and 
^//(O) ^/(O) = {2mV/h'^)'^j{0), the Andreev- and normal-reflection amplitudes a-f[j], 
for the spin- up (-down) quasiparticles are obtained as 

4n+Ai 



-■ --p- -_L 

(-1 - Ai - iz^[i]){-l -X2 + izm]) + (1 - Ai - iz^i]){-l + A2 - izm])n+n^(j)^(j)l 



(6) 



mi] 



-1 - A2 + izm]){l - Ai + iz^i]) + (-1 + A2 - izm]){-l - Ai + iz^[i])n+n-(j)-(j)* 
-1 - Ai - izm]){-l - A2 + iz^ii]) + (1 - Ai - iz^i]){-l + A2 - iz^^])n+n_(f)^(f)*^ 



(7) 



where zq 
given by 



mV 



220 

COS 9 s ' 



Ai 



cosg gT[i] 
cos 9s ks '' ^ 



S2^Amil_ The BCS coherence factors are 

cos 9s ks 



[i + JE^-\A^{e)\yE]/2, 



(8) 



\1 



E^-\Aii0)\yE]/2, 



(9) 



and n± = v±/u±. The internal phase coming from the energy gap is given by (f)± = 
[A±{0)/\A±{9)\], where A+{9) = A{9) (A_(6') = A(7r - 9)), is the pair potential experi- 
enced by the transmitted electronlike (holelike) quasiparticle. 

When 9 > sin~^(^) = 9ci total reflection occurs and the spin and charge current vanishes. 



In the space of 9, x in Fig. 0, the dotted line from the solution of the equation 9 = sin 
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/T+a:' 



where x = U/Ep, defines the boundary of the region (labeled as TR) where total reflection 



-l/£j_- 

superconductor, do propagate, the Andreev reflected quasiparticles, do not propagate. This 



occurs. When O^i > > sin (^) = 9^2 although the transmitted quasiparticles in the 



process is called virtual Andreev reflection (VAR) process |]T0|. In this case the spin and 
charge current do not vanish since a finite amplitude of the Andreev reflection still exists. 
For 6 < 6c2 Andreev reflection occurs. In Fig. ^ the solid line determined by the equation 



in the above inequality, i.e., 9 = sin~^ y irf > separates the region where the VAR process 
occurs (labeled as VAR) from the region where Andreev reflection occurs (labeled as AR). 
A symmetric branch that is not presented in the figure occurs for negative angles 6. 

According to the BTK formula the conductance for the charge current of the junction, 
ag^^^^{E,9), for up (down-)-spin quasiparticles, is expressed in terms of the probability am- 
plitudes am],b^i] as 



1 I -^2 I |2 IT |2 



Ai 

The tunneling conductance, normalized by that in the normal state is given by 



(10) 



a,{E) = a,^{E) + a,^{E), (11) 
^^Tm(^) = ^ r^l^^cos^a,^,^,(i?,^)PtU]gTU], (12) 

-TL/v J — 7r/2 



where 



= dO cos e[aN,{e)P^q^ + aN,{e)P^q^], (13) 

J —TV /I 
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where P|[|] = {Ep ± U)/2Ep is the polarization for up(down-)-spin. In the = limit the 
interface is regarded as a weak link, showing metallic behavior while for large zq values the 
interface becomes insulating. 

We consider the following cases 



a) In the {d^'^-yi + 2s)-wave case 

A{9) = Aicos[2{9- p)]+iA2, (15) 

where (5 is the angle between the normal to the interface and the x axis of the crystal. 

b) In the ((ia,2_y2 + idxy)-wsive case 

A(^) = Ai cos[2(^ - /?)] + iA2 sm[2{e - /?)], (16) 

where the angular form of the secondary component is obtained by the substitution of f3 in 

the {dx2-y2)-wave order parameter by /3 + 7r/4. 

III. TUNNELING CONDUCTANCE CHARACTERISTICS 

In Figs. and ^ we plot the tunneling conductance (yq{E) for different values of the 
exchange interaction x = U/Ep (a) Zq = 0, P = 0, (b) Zq = 2.5, (3 = 0, and (c) Zq = 2.5, 
j3 = tt/A. The pairing symmetry of the superconductor is d^^^y^ + is with Ai = Aq and 
A2 = O.3A0 in Fig. ^, dx2_y2 + id^y with Ai = Aq and A2 = O.SAq in Fig. For zq = 0, 
the subgap conductance is suppressed, with the increase of x, as in the case of a (ia.2_y2-wave 
superconductor [|To[] . 

In the {d^2_y2 + is)-wave case when the ferromagnet is normal metal (i.e., a; = 0), the 
boundary orientation is /5 7^ 0, and the barrier strength zq is large, a peak exists in the 
tunneling spectra in the amplitude of the secondary component due to the bound state 
formation. The peak height is maximum for (3 = n/A since the bound state is formed for all 
angles 9 and collapses to zero for /5 = 0. For the dj.2_y2 + id^y pairing state, for x = the 
tunneling conductance has residual values due to the formation of bound states. The bound 
state energies depend on the boundary orientation j3 as well as on the quasiparicle angle 6. 
The reduced height of the subgap conductance in the {dx2_y2 + idxy)-wauve case is explained 
from the discrete values of the angle 6 over which the bound state occurs as compared to 
the range of 9 values in the {dx2-y2 +is)-wave case [^. Also an enhancement appears in the 
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{d^2_y2 + idxy)-'wave state at x = 0, E = A^^^ for P = n/i due to the larger contribution to 
the spectra of the bound state at = 0{E = Ad^y).The same peak at A^^^ becomes more 
pronounced in a calculation including the self consistency ||16|| . 

In the Andreev-reflection process the incident electron and the Andreev reflected hole 
have wave vectors with opposite spins. In a normal metal the spin-up and spin-down wave 
vectors are equal and no spin effects occur in the Andreev reflection. However, in a ferro- 
magnet the wave vectors for spin-up and spin-down are different and this affects the Andreev 
reflection. In that case the Andreev reflected hole decays exponentially for large distance in 
the ferromagnet and there is no interference effect between electron and hole waves. More- 
over no pairs are transferred into the superconductor, and there is weak or no interference 
between the transmitted quasiparticles in the superconductor. In this sense the ferromagnet 
does not allow the quasiparticles to enter into the superconductor, and to experience the 
sign change of the pair potential, which is the main reason for the tunneling peaks. As a 
consequence the conductance peaks disappear when the exchange field gets very large. This 
is seen in Fig. ^(c), for the {dx2_y2 + is)-wave pairing state, where as the exchange field 
X increases the conductance peaks are reduced symmetrically. In the d^^^yi + id^y a dip 
develops within the subgap region as seen in Fig ^(c). The E = Q value is more sensitive 
to the exchange field (i.e., the Andreev reflection coefficient goes to zero faster) and the 
tunneling conductance for ii^ = 0, is suppressed more easily as the exchange field increases. 
In both pairing states the reduction of the subgap conductance is symmetric since the den- 
sity of states modulation within the subgap is not induced by spin-dependent effects, for 
example, a magnetic field. In that case we would expect an asymmetric evolution with the 
exchange field x since the effect of the magnetic field depends on the spin of the incident 
quasiparticle. This has been obtained in Ref. |T0| where the tunneling conductance in a 
ferromagnet/insulator/((i2.2_y2 + zs)-wave superconductor and also the effect of the mag- 
netic field if in a ferromagnet/insulator/((i2,2_j^2)-wave superconductor, is studied in order 
to identify the mechanism responsible for the splitting of the zero energy conductance peak 
in high-Tc superconductors. In the {dx2_y2 + 2s)-wave state for zq = 2.5, /5 = 0, and x = 0, 
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as seen in Fig. ^(b) there are no states within the subgap and <Jq{E) reduces to zero there. 
In the (i^2„y2 + id^y for zq = 2.5, /3 = as seen in Fig. ^(b) there are residual values within 
the subgap that are suppressed as x gets larger. For 2:0 = the evolution of the pairing state 
with the exchange field is similar in the two pairing states as seen in Figs. 0(a), §(a). 



IV. SUPPRESSION OF THE BOUND STATE ENERGIES 

We examine the evolution of the bound states with the exchange field x. The equation 
giving the energy peak level is written as 

(p^(l)ln+n_\E=Ep = ^-0- (17) 

When this condition occurs the Andreev- and normal-reflection amplitudes a^j^], feffj] for the 
spin up (down) quasiparticles are reduced to 

In Fig. ^(a) we plot the magnitude of the Andreev reflection amplitude for spin- up (-down) 
quasiparticle Re^|am]p as a function of the exchange field x for /5 = 7r/4, zq = 2.5, for the 
{dx2_y2 + is)-wave case. The corresponding magnitude of the normal-reflection amplitude 
is plotted in Fig. |^. The energy is equal to the amplitude of the s-wave component 
{E = O.3A0) for which bound states are formed for < ^ < 7r/2, when the ferromagnet 
is normal metal (i.e. x = 0). For 6 = 7r/4,7r/8, and x = 0, where Ai = A2, the Andreev 
reflection coefficient is equal to 1, and the normal-reflection coefficient is equal to 0, as 



obtained from Eqs. |T8| and |19[ In this case the conductance peak is due to the normal-state 
conductance in Eq. 14 that varies as 1/zq. As a; increases the amplitude of the Andreev- 



reflection decays to zero at a critical value Xc that depends from the angle 9 for which 
bound state occurs. The amplitude of the normal reflection increases with the exchange 
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field. The suppression of tfie Andreev reflection amplitude with x, explains the reduction of 
the conductance peaks as x increases. In the space of 6, x the critical exchange field Xc is 
defined from the separating line between the VAR region and the AR region in Fig. ^ For 
trajectories 6 that correspond to bound states the Andreev reflection vanishes within the 
VAR region. The critical exchange field Xc is maximum {xc = 1) when the bound state is 
at ^ = and is reduced to zero as 6 moves toward 6 = ±7r/2. This is also seen in Figs. |^(a) 
and H(b) where ioi 6 = ir/S, Xc = 0.75, while for 6 = it /4, Xc = 0.33. For a given value of the 
exchange field x, a bound state at 6 contributes to the cT(i?) only if x < Xc{0). This means 
that as X increases the range of bound states that contributes to the tunneling conductance, 
is reduced and the peaks are suppressed. 

In the {dx'2-y-2 + idxy)-v^a:Ve case when the ferromagnet is a normal metal, {x = 0) the 
bound states occur for discrete values of the quasiparticle angle 9, for fixed p. The Andreev 
reflection coefficient is equal to 1 for these values of 6. When the exchange field increases 
the Andreev (normal-)-reflection coefficient goes to 0(1) at a critical value Xc- This is seen 
in Fig. |^(a) and ^(b) for two different pairs of {E, 9) for which bound state is formed for 
X = 0, i.e., = 0,0 = 7r/4) and {E = O.3Ao,0 = 0) [§. The critical value of x for which 
the Andreev reflection coefficient goes to zero is independent from the pairing potential, and 
also from the energy of the bound state. It depends only on the angle 6. This is seen in Figs. 
^ and 1^, where for 6 = tt/4 the critical exchange field is Xc = 0.33 for both pairing states and 
for different values of the bound-state energy, i.e., E = 0.3Ao and E = correspondingly. 
The variation of the bound state angle 6 with Xc seen in Fig. ^ holds also for bound states in 
the {dx2_y2 + idxy)-^aMe state, and can be used to explain the suppression of the tunneling 
conductance at zero energy seen in Fig. ^c) with the exchange field x as follows. For 
E = 0.3Ao a bound state exists at = p[ for which Xc ~ 1 is maximum. For x < Xc the 
incident electrons are Andreev-reflected and the tunneling conductance has a finite value. 
The variation of the Andreev-reflection amplitude with x is seen in Fig. For the same 
energy another bound state exists at 6 = 7r/2, which does not contribute to the Andreev 
reflection since Xc is zero for this bound state. However, for E = one bound state is 
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formed for values of 9 close to — 7r/2 where Xc is close to zero and the Andreev reflection 
does not occurs. The other bound state at = occurs for = 7r/4 and contributes to the 
Andreev reflection up to x = 0.33, as seen in Fig. ^(a). For x > 0.33 the Andreev-reflection 
amplitude is zero and also the tunneling conductance is suppressed. This is seen in Fig. 
^(c) ioT E = and x = 0.4 (dotted line). Therefore the tunneling conductance at E = 
decays to zero more rapidly with x than the conductance at E = 0.3Ao. For a combination 
of {E, 9) for which no bound state is formed, the Andreev reflection amplitude is suppressed 
for all values of the exchange interaction x indicating that the exchange field mainly affects 
the bound states. 

V. MAGNETIC FIELD EFFECTS 

In this section we describe the effect of the external magnetic field H in the spectra for 
different values of the exchange field x. We will see that since the effect of the magnetic 
field depends on the spin, the evolution of the tunneling spectra with x is asymmetric. The 
tunneling conductance is given by 

a,{E) = a,^{E - ^ibH) + a^^ {E + ^ibH). (20) 

In Fig. 0(a) and |^(b) the tunneling conductance crq{E) is plotted for fixed magnetic field 
(ibH/Aq = 0.2, and barrier strength zq = 2.5, for different values of the exchange interaction 
X. The pairing symmetry of the superconductor is dx2-y2 +is and d^i-yi +idxy, respectively. 
The orientation of the superconductor is chosen as /3 = tt/A. 

In the absence of the exchange interaction [x = 0) the magnetic field splits symmetrically 
the tunneling spectrum that is a linear superposition of the spectra for spin up (down) 
quasiparticles. The amplitude of the spliting depends linearly on the magnetic field H. For 
the case of fisH/Ao = 0.2, seen in Fig. |^ the spin- up (-down) part of the spectra partially 
overlap while for larger values of the magnetic field the spin-up and -down branches are well 
separated. In the latter case the left(right) branch of the spectra originates from spin-up(- 
down) quasiparticle spectra <Jqy{E — jj,BH)[aq^{E + hbH)]. 
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For the dx2_y2 + is-wave case the condition for the formation of bound states is shghtly 
modified under the presence of magnetic field to \E — ^bH\ = A2, for the spin-up region, 
and \E + ^bH\ = A2, for the spin-down, from the corresponding \E\ = A2 in the absence of 
any field. So the multiplication of the {d^^-yi + zs)-wave pairing state and the presence of 
magnetic field results into the appearance of four peaks in the conductance spectra, which 
in the limit of x = have equal heights. 

The main effect of the polarization is the imbalance in the peak heights for E positive 
and negative. The ratio of the peaks for positive and negative energy is proportional to 
the exchange field of the material. This can be extracted from the different evolution of 
the Andreev- and normal-reflection coefficients for spin-up and -down quasiparticles with 
the exchange field seen in Figs. § and ^ Note that although the bound state energies are 
modified in the presence of the magnetic field, the analysis concerning the above Figures 
still holds for the modified energies. For a given energy and angle 9 for which bound state 
occurs the quantities aq^{E, 9), a^^^E, 9), have different values causing the asymmetricity in 
the peak heights for the spin-up and -down part of the spectrum. The asymmetricity in the 
Andreev reflection coefficient can also be seen in Fig. for spin up (down) quasiparticles 
as solid (dotted) line, as a function of the energy E/Aq for fixed exchange field x = 0.6, 
and the pairing symmetry of the superconductor is d^2_y2 + is. The same characteristic 
appears in the normal reflection coefficient that is plotted in Fig. H(b). The peaks in the 
Andreev-reflection coefficient are formed at the bound state energies and due to the finite 
exchange interaction are suppressed from the unit. The same result is plotted in Fig. ||(a) 
and §(b), for the ((i^2_y2 -|- idxy)-"wave pairing symmetry. Here the bound state is formed 
for a particular value of 6* = 0. Other reasons for the asymmetricity of the spectra for the 
spin up(down) quasiparticle are the factors P|[|] and that appear in the definition of the 
tunneling conductance Eq. (|T2|). 
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VI. C-AXIS TUNNELING 



In the preceding sections we discussed the tunnehng effect in two dimensional models. 
In this section we discuss the tunneling effect along the c axis that takes into account three- 
dimensional effects. A semi-infinite double layer structure is assumed and the volume of the 
integration is taken as the three dimensional half sphere. The interface is perpendicular to 
the z axis and is located at z = as seen in Fig. Suppose that an electron is injected 
from the ferromagnet with polar angle Q and azimuthal angle 0. The electron like (hole) like 
quasiparticle will experience different pair potentials A^^/(^^+) [A^^'(^^_)], where 6^ = 6 and 
9- =71 — 6, and the quantities p, p denote spin indices. The coefficients of the Andreev and 
normal reflection are obtained by solving the BdG equations under the following boundary 
conditions 

vl/(r)|,=o_ = *(r)|.=o+ (21) 

-^U=0- - ^1-0, - ^^(r)U=o_ (22) 

Using the obtained coefficients the tunneling conductance is calculated using the formula 
given in the preceding sections, 

I i-n/2 i-2n 

^9Tm(^) = ^/ / <^ose sin eaq^^^^{E,e,(j))Pm]qm]ded(j), (23) 

where 

Rn= smecose[aN^.{e)P^q^ + (7N^{e)Piqi]ded(j), (24) 





4Ai 



2 + ^2 



(25) 



The pairing potentials are given by 
a) In the {d^^-yi + is)-wave case 



A(e,0) = Aicos20 + iA2, (26) 
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b) In the {dx2_y2 + idxy)-'wa:Ve case 

A(6',</)) = Aicos2(/) + 2A2sin2(/). (27) 

Figures |ll](a) and |ll](b) show the calculated conductance spectra for various exchange po- 
tentials. Unlike the case where the interface is perpendicular to the x axis the tunneling 
spectra is similar for the dx2_y2+is, dx2-y2+idxy cases. Also a subgap region is formed within 
the energy gap due to the nodeless form of the order parameter. The conductance peaks are 
absent since the transmitted quasiparticles do not feel a sign change of the pair potential 
on the Fermi surface. The case of c-axis tunneling has been treated experimentally using 
scanning tunneling microscopy or point-contact spectroscopy, for the case of normal metal 



/ superconductor junction ||T^|18|. No zero energy peak and a clear V^-like line shape for the 
c-axis tunneling conductance has been observed. The calculated spectra for the dx^_y2 + is, 
dx'2-y2 + idxy states is flattened out for E < A2 that contradicts the experimental data. 
The ((ia,2_y2)-wave order parameter does not change sign and is not suppressed at the c-axis 
surface. So the attractive interaction in the subdominant pairing channel is small relative to 
the dominant and a transition to a state breaking T does not happen. On the other hand a 
mixed order parameter such as the dx2_y2 + s may exists due to the orthorhombic distortion 
of the lattice, as seen in the c-axis Josephson experiments. 



VII. CONCLUSIONS 

We calculated the tunneling conductance in ferromagnet / insulator / superconductor, 
junction using the BTK formalism. We assumed two possible pairing potentials for the 
superconductor that break the time-reversal symmetry, i.e., dx2_y2 + is, dx2-y2 + id^y. The 
evolution of the spectra with the exchange field is the same for zq = {] but different in 
the tunneling limit where zq is large, and can be considered as a probe for time reversal 
symmetry broken pairing states. The weak Andreev reflection within the ferromagnet results 
in the suppression of the tunneling conductance and eliminates the resonances due to the 
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anisotropy of the pair potential. The evolution of the tunneling conductance within the gap 
is symmetric since the splitting is not induced from spin- dependent effects, for example, a 
magnetic field, but from the nodeless form of the pairing potential. 

We also derived that the condition for a bound state at angle 9 to contribute to the 
Andreev-reflection and hence to the tunneling conductance, for a given value of the exchange 
interaction x is x < Xc- Xc is the critical exchange field for which the Andreev refletion 
coefficient goes to zero and is given from the separating line between the VAR region and 
the AR region in the space of 6,x. This condition was used to explain the suppression of 
the conductance around E = with the exchange field in the ((ij.2_y2 + idxy)-waye state. 

The magnetic field splits linearly the tunneling spectra, and the exchange potential 
induces an imbalance in the peak heights for positive and negative energies. The asym- 
metricity in the peak heights originates from the different evolution of the Andreev- and 
normal-reflection amplitudes at the bound state energies with the exchange field. 

The c-axis tunneling from ferromagnet to superconductor does not show any differences 
between T-broken pairing states since the transmitted quasiparticles experience the same 
sign of the pairing potential. 

Throughout this paper the order parameter is not calculated self-consistent ly. However, 
since the characteristics of the tunneling conductance depend mainly from the angular part 
of the pairing potential, the essential results are expected to change only quantitatively when 
the suppression of the order parameter near the surface is taken into account. 
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electron-like 
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hole-like 
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FIG. 1. The geometry of the ferromagnet / insulator / singlet superconductor interface. The 
vertical line along the y axis represents the insulator. The arrows illustrate the transmission and 
reflection processes at the interface. 6 is the angle of the incident electron and the normal, 9a_ is 
the angle of the reflected-hole trajectory and the normal, and 9s is the angle of the transmitted 
quasiparticle and the normal. Note that 9 is not equal to 9a since the retroreflection of the Andreev 
process is lost. In the Andreev-reflection process an electron with spin up is Andreev reflected as 
a hole with spin down and normally reflected as an electron with spin-up. 
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FIG. 2. The regions of the space of 9, x TR, where total reflection occurs, VAR where Andreev 
reflected quasiparticles do not propagate, while the transmitted quasiparticles propagate, AR where 
Andreev reflection occurs. 
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FIG. 3. Normalized tunneling conductance crq{E) as a function of E/Aq for x = (solid 
line), X = 0.4 (dotted line), x = 0.8 (dashed line), and x = 0.999 (long-dashed line), for different 
orientations (a) zq = 0, (3 = 0, {h)zo = 2.5, (3 = 0, (c) zq = 2.5, /? = 7r/4. The pairing symmetry of 
the superconductor is d^2_y2 + is with Ai = Aq, A2 = O.SAq. 
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FIG. 5. a) The magnitude of the Andreev reflection coefficient Re^|ai-[j]|^ as a function of the 
exchange field x, for spin- up (-down) quasiparticles, solid (dotted) line, for ^ = 7r/4 and 6 = tt/8, at 
the bound-state energy E = O.SAq. The pairing symmetry of the superconductor is d^.2_j^2 +is and 
/3 = 7r/4. b) The corresponding magnitude of the normal-reflection coefficient |fe-f[j]P as a function 
of the exchange field x. 
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FIG. 6. The magnitude of the Andreev reflection coeflicient Re^|a|[j]p (a) and nor- 
mal-reflection coefficient (b) as a function of the exchange fleld x, for spin-up (-down) 
quasiparticles, sohd (dotted) hne, for two pairs of {E,6), for which bound states occur i.e. 
{E = 0,9 = '/r/4), and {E = O.SAq,^ = 0) in a superconductor with {d^2_y2 + idxy)-wave pairing 
and /3 = 7r/4. 
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FIG. 7. Normalized tunneling conductance CTq{E) as a function of E//S.Q for = (solid 
line), X = 0.2 (dotted line), x = 0.4 (dashed line), and x = 0.999 (long-dashed line), for zq = 2.5, 
/3 = 7r/4, in the presence of an external magnetic field hbH/Aq = 0.2. The pairing symmetry 
of the superconductor is (a) dr^2_y2 + is with Ai = Aq and A2 = O.SAq. (b) dr^2_y2 + idxy with 
Ai = Ao and A2 = O.SAq. 
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FIG. 8. a) The magnitude of the Andreev reflection coefficient Re^|a-|-^]p as a function of 
E/Aq, for spin- up (-down) quasiparticles, solid (dotted) line, for 6 = vr/8, for exchange field x = 0.6, 
fisH /Aq = 0.2. The pairing symmetry of the superconductor is d^2_y2 + is and (3 = n/A. b) The 
corresponding magnitude of the normal-reflection coefficient Ifeff^jP as a function of E/Aq. 
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FIG. 9. a) The magnitude of the Andreev reflection coefficient Re^|a-|-[j]p as a function of 
E/Aq, for spin-up (-down) quasiparticles, solid (dotted) line, for ^ = 0, for exchange field x = 0.6, 
fisH/Ao = 0.2. The pairing symmetry of the superconductor is dj.2_y2 + idxy, and (3 = 7r/4. b) 
The corresponding magnitude of the normal-reflection coefficient as a function of E/Aq. 
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FIG. 10. The figure illustrate the transmission and reflection processes of the quasiparticle at 



the interface of the junction with xy plane interface. 
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FIG. 11. Normalized tunneling conductance CFq{E) as a function of E/Aq for x = (solid 
line), X = 0.4 (dotted line), x = 0.8 (dashed line), and x = 0.999 (long dashed line), for zq = 2.5. 
The interface is perpendicular to the z axis. The pairing symmetry of the superconductor is (a) 
dx'i-yi + is with Ai = Aq and A2 = O.SAq. (b) d^2_y2 + idxy with Ai = Aq and A2 = O.SAq. 
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